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Abstract. Near-infrared observations of line emission from excited H2 and in the continuum are reported in the 
direction of the Orion molecular cloud OMC1, using the European Southern Observatory Very Large Telescope 
UT4, equipped with the NAOS adaptive optics system on the CONICA infrared array camera. Spatial resolution 
has been achieved at close to the diffraction limit of the telescope (0.08" - 0.12" ) and images show a wealth of 
morphological detail. Structure is not fractal but shows two preferred scale sizes of 2.4" (1100 AU) and 1.2" (540 
AU), where the larger scale may be associated with star formation. 
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1. Introduction 

Near infrared (NIR) observations, performed with the 
European Southern Observatory Very Large Telescope 
(ESO-VLT), are described which provide high spatial res- 
olution images of the Orion Molecular Cloud (OMC1; dis- 
tance 460 pc (Bally et al. 2000)). Data are presented in the 
v=l-0 S(l) line of H2 at 2.121/xm and in continuum emis- 
sion at nearby wavelengths. Earlier NIR imaging of this 
zone may be found for example in Allen&Burton 1993; 
Brand 1995; Schild et al. 1997; McCaughrean&MacLow 
1997; Chrysostomou et al. 1997; Chen et al. 1998; Stolovy 
et al. 1998; Schultz et al. 1998; Lee&Burton 1999; Tedds 
et al.1999; Salas et al. 1999; Vannier et al. 2001 (V2001); 
Gustafsson et al. 2003 (G2003); Kristensen et al.2003 
(K2003). The nature of the Orion Nebula Cluster, in 
which OMC1 is embedded, is reviewed in O'Dell 2001 and 
Ferland 2001. As discussed there, OMC1 consists of dense 
fragments of gas and dust, in part situated within the HII 
region, which itself constitutes the visible Orion Nebula. 
The interest in OMC1 and the surrounding zone stems 
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from widespread ongoing star-formation in this region 
(e.g. Luhman et al. 2000), exemplified through the pres- 
ence of many protostars, outflows and larger scale flows, 
for example the explosive Becklin-Neugebauer-IRc2 com- 
plex (BN-IRc2) in OMC1 (see O'Dell 2001; Allen&Burton 
1993; Doi et al. 2002; O'Dell&Doi 2003). 

In the present work we use the VLT-UT4 with NACO, 
an association of the new adaptive optics system NAOS 
(Nasmyth Adaptive Optics System) and the CONICA 
(COude Near IR Camera) infrared array camera (Lenzen 
et al. 1998; Rousset et al. 2000; Brandner et al. 2002). This 
yields, in part of the full image, a spatial resolution two 
to three times higher than that achieved previously with 
the HST (Stolovy et al. 1998; Schultz et al. 1998) or with 
4m-class telescopes using adaptive optics e.g. the Canada- 
France-Hawaii Telescope (V2001; G 2003; K2003). 

Our present data allow us to address the issue of the 
scale sizes of clumps of excited H 2 in OMC1 with signif- 
icantly greater precision and over a region within OMC1 
roughly an order of magnitude larger in the plane of 
the sky than in earlier work (V2001). In this connection, 
V2001 showed that H2 emitting regions did not exhibit 
fractal structure but showed a preferred scale of between 
3 and 4 10~ 3 pc. Here we corroborate this finding and 
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Fig. 1. H2 emission in the 1-0 S(l) line at 2.121/xm, with continuum emission subtracted. Image size: 44.5"x26". The 
cross indicates the position of the star TCC016 (05 h 35 m 14?91, -05°22 , 39'/31 (J2000)) and a,b,c identify objects shown 
in more detail in subsequents. The SE and ESE frames, referred to in the text, are respectively the regions to the west 
and to the east of TCC0016. 



also identify a second preferred scale size of about half 
this value, in a neighbouring region of OMC1, to the west 
of that studied in V2001. 

2. Observations and data reduction 

Observations of OMC1 were performed with the VLT UT4 
during guaranteed NACO time on the 21 st of November 
2002, for about 4.5 hours in total. The seeing was vari- 
able during this night ranging from 0.65" to 1.13". Two 
regions have been observed located south-east (SE frame) 
and east-south-east (ESE frame) from BN and IRc2, with 
the star TCC016 (Fig. 1) common to both frames. The IR 
wavefront sensor (WFS) was used for the SE frame using 
BN (111^—8) as the adaptive optics (AO) reference star. 
The average seeing was —0.85" and the resulting FWHM 
of the PSF measured on several stars in the SE field was 
—0.08". The ESE frame was recorded using the visible 
WFS locked on the star TCC016 (m y -14) which, com- 
bined with an average seeing of —1.03" , led to a measured 
FWHM of the PSF of -0.12" . The Strehl ratio was clearly 
higher in the SE frame, giving nearly full correction by 
the AO system to the diffraction limit of —0.06". Using 
camera mode S27, the field of view was 27.6" x27.6" and 
the pixel scale was 0.027", a size sufficient to satisfy the 
Nyquist sampling criterion. 0.027" corresponds to 12 AU 
at the distance of Orion. 

Filters were used at 2.121±0.011/mi (NB212), which 
contains the H 2 S(l) v=l-0 line, IB224, at 2.24±0.03/zm 
and IB227 at 2.27±0.03/itm. IB224 includes the H 2 1-0 S(0) 



line at 2.223/xm and the H 2 2-1 S(l) line at 2.247/xm, as 
well as [Fell] lines at 2.218 and 2.242^m. IB227 includes 
the H 2 2-1 S(l) line and the [Fell] line at 2.242^m. The 
AutoJitter mode was used, that is, at each exposure, the 
telescope moves according to a random pattern in a 6" x 6" 
box. Cross-correlation was used to recenter the images at 
«0.15 pixel. For the ESE frame, an empty background 
sky has been recorded for the purpose of sky subtraction. 
For the other less crowded SE frame, the background sky 
has been created by performing a median filtering of the 
set of 32 randomly jittered individual frames. The total 
exposure time on object for each filter was 800 seconds 
for the SE frame and 1500 seconds for the ESE frame, 
with exposures of 10 seconds for each data acquisition. 

To obtain an image in the H 2 v=l-0 S(l) line, the 
IB244 continuum image was subtracted from the NB212 
image after division by a factor of 2.5 as derived from 
measuring fluxes of a dozen of stars in both filters. The re- 
sulting image is shown in Fig^ The contribution of v=T 
S(0) and v=2-l S(l) in IB224, known to be small from 
spectroscopy, was ignored in image processing. By sub- 
traction of the image in the IB227 filter from that in the 
IB224, using a conversion factor obtained as before from 
stellar fluxes, it was possible to obtain an image in H 2 S(0) 
v=l-0, less the weaker emission in v=2-l S(l). The image 
in v=l-0 S(0) shows essentially the same structure as that 
in Fig^ but is 4 to 5 times weaker. An example of data 
in v=l-0 S(0) is shown in Fig|3J:. 

Strong continuum emission is seen for example in the 
northern part of the ESE region and is shown in Fig|5J 
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Fig. 2. Continuum emission at 2.27±0.03/mi in the ESE 
field. Image size: 27.6"x23.4" . The position of TCC0016 is 
marked with a cross. 
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Fig. 3. a). A detailed view (3.5"x5.0") of object a (see 
FigEJ in v=l-0 S(l) H 2 emission at 2.121/xm. b). of object 
b (2.5"x5.0") in the same line. c). Detail (2.5"x5.0") of 
object b in the 1-0 S(0) H2 emission line at 2.22/im. 



which displays emission in the IB227 filter. Data in Fig|21 
also include weak emission from v=2-l S(l) and [Fe II] 
(Schultz et al.1999). Continuum emission is also seen at 
this position in the data of Schultz et al. 1999. 

3. Discussion of Results 

3.1. Origin of the NIR emission 

H2 emission in OMC1 arises from heating in J-type and C- 
type (magnetic) shocks (e.g. Draine et al. 1983, Pineau des 
Forets et al.1988, Smith&Brand 1990, Kaufman&Neufeld 
1996a,b, Timmerman 1998, Wilgenbus et al. 2000, V2001, 
Le Bourlot et al. 2002, K2003, G2003) and from 
photon excitation in photodissociation regions (PDRs) 
e.g. Storzer&Hollenbach 1999, Sternberg&Dalgarno 1989, 
Black&vanDishoeck 1987, Black&Dalgarno 1976. 

Shocks in the region observed arise from one or more 
large-scale outflows, thought to originate from the BN- 
IRc2 region (O'Dell&Doi 2003, Doi et al. 2002, O'Dell 




Fig. 4. A detailed view of object c (see Fig^l in v=l-0 
S(l) H2 emission at 2.121/xm. Image size: 7.5"x5.0". 



2001 and references therein). Shocks in the zone observed 
may also arise from local outflows associated with proto- 
stars buried within emitting clumps of gas (G2003). 

PDRs are generated through the action of the 
Trapezium stars, to the south of the present region, of 
which the dominant contributor is ^Ori C. This star gen- 
erates far-UV radiation fields of >10 5 times the standard 
interstellar field over much of OMC1, and is in addition the 
chief contributor to the HII region in which dense clumps 
in OMC1 may be bathed. The BN object is a young and 
massive B-star (Gezari et al. 1998) and may also provide 
a source for PDR excitation. 

The relative importance of PDR and shock excitation 
in the strongest emitting regions has been discussed in 
detail in K2003. The conclusion is that in the subset of 
clumps studied in that work, which are those to the east 
of TCC016 in the ESE image, C-type shocks are the major 
contributor to H2 excitation. However, towards the fringes 
of bright clumps, J-type shock and PDR excitation by 
6' 1 Ori C take the place of C-type shock excitation, with 
shock and PDR excitation contributing roughly equally. 
In all cases of very bright emission, both in the cores of 
clumps and at their fringes, the densities are high, exceed- 
ing several times 10 7 to >10 8 cm~ 3 . Excited gas, formed 
in numerous, individually unresolved shocks, rapidly accu- 
mulates into cold compressed zones and excited H2 is the 
very rapid progenitor of cold dense gas. Thus the struc- 
ture of excited H2, discussed in section 3.2 below, provides 
a measure of the structure of cold H2 clumps, where the 
latter is the gas from which stars will ultimately form. 

In other regions, which are presumably less dense, 
the correspondingly larger shock structure appears to be 
clearly resolved in H2 emission. This is illustrated by data 
in Fig. |31 El and El Numerous objects in these figures re- 
semble bow shocks. The general appearance of this figure 
suggests that the region may be permeated by supersonic 
turbulence giving rise to shocks in all directions. The emis- 
sion may also in part be due to complex density structure 
with PDR excitation. At all events, our NACO data allow 
us to specify the width of the filamentary regions con- 
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taining excited gas. For example in Fig|SJ showing the 
western part of object c (Fig0J and of object a (Fig|3J) in 
greater detail, widths of some filaments appear marginally 
resolved. If these scales of tv40-50 AU are shock widths, 
these data provide important constraints on shock models 
(Le Bourlot et al. 2002, Wilgenbus et al. 2000), in which 
parameters of density, magnetic field and shock speed de- 
termine the shock width. If the features are PDRs, the 
widths of structures are again valuable parameters, pro- 
viding a good indicator of the gas density (K2003, Lemaire 
et al. 1996, Sternberg&Dalgarno 1989). An additional con- 
straint to shock or PDR models is provided by data for the 
v=l-0 S(0) line, see Fig|!J;. In future work, we will pur- 
sue modeling of these regions, including constraints arising 
from absolute brightness in the H 2 emission lines. 

The origin of the diffuse continuum emission (plus 
[Fell] and some weak H2 emission - see above), for ex- 
ample in Fig|2 is most likely reflected light from dust. 
In this connection, there is little shock excitation where 
strong continuum emission is observed, as a comparison of 
Fig and [21 reveals. Thus little mechanical energy is being 
injected into this zone. Emission could arise from far-UV 
photon heating of very small dust particles, as for exam- 
ple in NGC7023 (Lemaire et al. 1996), with 6 1 Oii C and 
BN as sources in the present case. There is however strong 
circular polarization in the K-band continuum in the zone 
shown in Fig|21 (Chrysostomou et al. 2000). Those authors 
ascribed the origin of circular polarization to scattering 
from oblate grains oriented in magnetic fields. There is 
in fact striking spatial correlation, lying within the ob- 
servational resolution of the circular polarization data, 
between the regions of strong circular polarization re- 
ported in Chrysostomou et al. 2000 and the continuum 
emission seen in Fig |3] and at other positions throughout 
the SE and ESE images. If we include the linear polar- 
ization data of Geng 1993, then effectively in all regions 
in which continuum emission is observed, cither circular 
or linear polarized K-band radiation may be found. This 
lends strong support to the hypothesis that the observed 
diffuse emission in the continuum is due to scattering from 
dust. According to Chrysostomou et al. 2000, the source 
of illumination is not the Trapezium stars, but is intrinsic 
to OMC1 through some unidentified source. 

Weak [Fell] emission is present for example in the re- 
gion shown in Fig|2 (Schultz et al.1999) and most likely 
arises from photoionization and excitation by 8 1 Oii C in 
the HII plasma as suggested in Schultz et al. 1999. In this 
model, [Fell] emission is projected over the K-band con- 
tinuum emission discussed above. 

3.2. Analysis of small scale structures 

V2001 gave the first clear indication that the size distribu- 
tion of hot H2 clumps in Orion was not fractal, using data 
for a small region of OMC1, 12.8" x 12.8", in the south- 
ern part of the ESE region of the present work. In or- 
der to analyze the global nature of small scale structure, 



V2001 used area-perimeter, Fourier, brightness histogram 
(Blitz&Williams 1997) and direct measurement analyses. 
All techniques revealed that structure observed in v=l- 
S(l) H2 emission was not fractal and all yielded the 
same preferred scale size. Here we choose the brightness 
histogram method for its simplicity and ease of interpre- 
tation. 

In brief, this method involves counting in the original 
image the number of pixels N(B n ), binned into some small 
normalized brightness range, as a function of the pixel 
normalized brightness, B„. B„ is defined as the value of 
brightness in any pixel divided by the value in the bright- 
est pixel in the image. The resolution of the image is then 
degraded by boxcar averaging over some chosen number 
of pixels on a square and a new set of N(B„) and B„ cal- 
culated. This process is repeated, successively degrading 
the resolution of the image. Plots of N(B n )/NTot vs. B„, 
where Nrot is the total number of pixels, yield a set of 
curves which remain unchanged for an image with no pre- 
ferred scale, that is, a fractal image, as the resolution of 
the image is degraded. By contrast, for an image with a 
preferred scale, such plots should change in form as the 
resolution is degraded to the point at which any preferred 
scale has been washed out. Beyond this critical smoothing, 
the form of such plots should remain constant, mimicking 
a fractal. The amplitude of the critical smoothing is a di- 
rect measure of the preferred range of scale in the image. 
In all cases stars were masked out in the images before the 
brightness histogram analysis was performed. 

A brightness histogram analysis for the ESE image in 
H2 is shown in FigEl We estimate that the preferred scale 
for regions of excited H 2 is encountered when averaging is 
performed over 90±5x90±5 pixels, that is, at a scale of 
2.4±0.14" (1075 AU). An exactly similar analysis on the 
SE region, that is, to the west of TCC016 in Fig^ yields 
a preferred scale of 1.2±0.14" (540 AU). If we choose a re- 
gion which samples approximately equally both the ESE 
and SE regions and perform a brightness histogram anal- 
ysis, we encounter two scales present together, of 2.4" and 
1.2". We conclude therefore that the ESE region contains 
clumps of material of about twice the scale of the ma- 
terial in the SE region, with a well-defined demarcation 
between them. The same analysis has been applied to the 
image in the continuum, shown in Fig[2 This yields a 
scale of 2.15±0.28". This indicates that dust, and by im- 
plication, cold unexcited gas in this region is clumped on 
the same scale as in the shocked southern region in ESE. 
The scale size, if one is present, of continuum emission 
in the SE region could not be determined using the tech- 
nique described here, but is larger than 2.7 arcseconds. 
Data covering a more extensive field in OMC1 may yield 
information in the future. 

The scale-size identified in V2001, for a subset of the 
ESE region, was 2.0±0.1", in substantial agreement with 
the present work. As described in detail in V2001, the 
brightness of emission of H2 may be used in conjunction 
with C-type magnetic shock models to show that in bright 
emissions regions the dense gas, swept up by the shock, 
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Fig. 5. a). Detail (1.56"x2.40" ) of the NW part of FigH 
and b). of Fig|3t (2.24"x2.40" ) showing narrow filamen- 
tary structure. 

will achieve densities of several times 10 7 cm -3 , as the 
post-shock gas cools to 10K. Related data for the v=2-l 
S(l) emission line, discussed in K2003, show that densities 
may be locally even higher, exceeding 10 8 cm" 3 in some 
regions. Our present estimates of clump size add weight 
to the conclusion in V2001, based upon Jeans length con- 
siderations, that sufficient gas has accumulated through 
shock induced compression to form gravitationally unsta- 
ble clumps which are potential sites of star formation, 
in the ESE region. This shock-induced mechanism would 
not appear to be presently effective in hastening star- 
formation in the SE region, since we find that the clump 
masses are typically an order of magnitude smaller than 
in ESE and comprise no more than 0.01 to 0.02 M Q if we 
assume densities of < 10 8 cm~ 3 given the observed size of 
1.2". In connection with the interpretation of our data, 
results in G2003 appear to show that protostars may al- 
ready have formed in certain of the zones very bright in H2 
emission. In some regions we may therefore be observing 
the aftermath of very recent star formation rather than 
the period just prior to gravitational collapse. 
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